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Abstract. Some active stars show a so-called flip-flop phenomenon in which the main spot activity periodically 
switches between two active longitudes that are 180° apart. In this paper we study the flip-flop phenomenon by 
converting results from dynamo calculations into long-term synthetic photometric observations, which are then 
compared to the real stellar observations. We show that similar activity patterns as obtained from flip-flop dynamo 
calculations, can also be seen in the observations. The long-term light-curve behaviour seen in the synthesised 
data can be used for finding new stars exhibiting the flip-flop phenomenon. 
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1. Introduction 

In many active stars the spots concentrate on two perma- 
nent active longitudes which are 180° apart. In some of 
these stars the dominant part of the spot activity changes 
the longitude every few years. This so-called flip-flop phe- 
nomenon was first reported by Jetsu et al. ifTMTl [T^Mjl 
in the single, late type giant FK Com. Berdyugina & 
Tuominen l|1998|l reported periodic flip-flops between per- 
manent active longitudes in four RS CVn binaries. Their 
results were confirmed in the case of II Peg by Rodono 
et al. ((5000). The persistent active longitude structures 
and flipping between two active longitudes have also been 
reported over the years based on photometric observa- 
tion (e.g. Berdyugina et al. Korhonen et al. I^OH^ 
Jarvinen et al l2005|l and on Doppler images (Berdyugina 
et al. 119981 Korhonen et al. I2001|l . After its discovery in 
cool stars, the flip-flop phenomenon has also been reported 
in the Sun (Berdyugina & Usoskin|2003)). A review on the 
flip-flop phenomenon in cool stars and the Sun is given by 
Berdyugina H2004|l . 

In order to explain this phenomenon, a non- 
axisymmetric dynamo component, giving rise to two per- 
manent active longitudes 180° apart, is needed together 
with an oscillating axisymmetric magnetic field. Fluri & 
Berdyugina (200J| suggest also another possibility with a 
combination of stationary axisymmetric and varying non- 
axisymmetric components. Unfortunately, no simple dy- 
namo mechanism is yet known for such a configuration. 
There are models with anisotropic a-efFect or a weak dif- 
ferential rotation, which could produce non-axisymmetric 
components, but only recently Moss £2004. .2005,1 reported 
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coexisting mixed components with a differential rotation 
that depends on the distance to the rotation axis. This 
differential rotation configuration is a very plausible state 
for fast rotators. Also, a weak non-axisymmetric field co- 
existing with a dominant axisymmetric field, assuming a 
solar-like rotation law, was found by Moss (1999). This 
solution was not analysed for flip-flops, as the possibility 
of them being present on the Sun was not being discussed 
at that time. Flip-flop solutions for a rotation law similar 
to the solar one and anisotropic a have also been found 
by Elstner & Korhonen 

According to the calculations by Elstner & Korhonen 
1)2005(1 the shift of the spots in a flip-flop event is 180° only 
in some cases, mainly for the stars with thin convective 
zones. In stars with thick convective zones they found a 
shift that is close to 90°. Similar results were reported by 
Moss ("2005) . Recently, Olah et al. f2005') re-analysed some 
of the old photometric observations of FK Com and found 
a flip-flop event in which spots on both active longitudes 
vanished briefly, and one of the new spots appeared on an 
old active longitude and the other one 90° away from that. 
This is the first evidence suggesting that flip-flops where 
the spots shift only by 90°can also occur. 

In this paper we use the model calculations presented 
by Elstner & Korhonen 1(2005(1 and convert them into syn- 
thetic photometric observations. This is used to investi- 
gate the expected long-term photometric behaviour of ac- 
tive stars showing the flip-flop phenomenon. Hopefully, 
this will help us in identifying new stars exhibiting this 
intriguing phenomenon. At the moment only few stars 
showing it are known and no statistically significant cor- 
relation between the stellar parameters and the flip-flop 
phenomenon can therefore be carried out. 
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2. Model 

The model consists of a turbulent fluid in a spherical shell 
of inner radius rin and outer radius rout • 
We solve the induction equation 



d{B) 
dt 



curl(ag o [B) — ?7TCurl(B)), 



(1) 



in spherical coordinates (r, 0, Lp) for an a^f2-dynamo. Here 
we used a solar-like rotation law in the corotating frame 
of the core 



1 + erf 



(2) 



where i7s = f^eg — acos^O with fleq/^ir — 460.7nHz, 
J]c/27r = 432.8nHz and a/27r = 125.82nHz are used. 
Normalising length with stellar radius i?* and time with 
diffusion time td = Rl/i] we can define the dimensionless 
dynamo numbers 



and 

Cn = i^eg - n,) ■ Rl/r,. 



(3) 



(4) 



For Ri, = i?0, the diffusivity 77 = 5 ■ 10^^cm-^s~^ and JIq — 
1 we get Co = 172. Notice that a definition with Q^q — 
^poi gives Co — 780. Because of these different definitions 
of Cn in the literature we describe the strength of the 
differential rotation in our models with the value of ^Iq. 

In order to identify the lifetime and the maximal possi- 
ble pole to equator difference of the angular velocity for a 
flip-flop solution also for models with isotropic a, we per- 
formed several calculations with models similar to those 
presented in Moss (,,2005.) . Here we used the differential 
rotation law (Eq. |2J) with fig — —ar^sin^d and Qc — 0. 
The same normalisation was used. 

Only the symmetric part 

arr = aocos6'(l. — 2cos^^^) 
o:ee = aocos0(l. — 2sm^9) 
a.^^ = aoCOsO 
are — cegr = 2aocos^9sm9 



(5) 



of the a-tensor is included (cf . Riidiger et al 120021) . For 
the isotropic model are — oter = and the diagonal terms 
Oiii = a^cosd. In order to saturate the dynamo we choose 
a local quenching of 

(6) 
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The inner boundary is a perfect conductor at rin — 0.3 
and the outer boundary resembles a vacuum condition, by 
including an outer region up to 1.2 stellar radii into the 
computational grid with 10 times higher diffusivity. At the 
very outer part the pseudo vacuum condition (tangential 
component of the magnetic field and vertical component of 
the electric field vanish on the surface) is used. In order to 
see the influence of the thickness of the convection zone 



Table 1. Critical Cq for the symmetric and antisym- 
metric components with azimuthal wavenumbers m=0,l. 
The second number gives the oscillation period (for m=0) 
or the migration period (for m=l) in diffusion times (30 
years) . 



name 


SO 


SI 


AO 


Al 


thick 


9.7 0.24 


9.0 0.15 


9.7 0.24 


9.0 0.15 


thin 


11.9 0.18 


10.2 0.27 


11.9 0.18 


10.2 0.27 



we have chosen = 0.4 for a thick (results shown in 
Fig-Et) and Tc = 0.7 for a thin (Fig.lSJj) convection zone. 
Below the convection zone from rin to Tc the diffusivity 
was reduced by a factor 1000 and a was set to zero. The 
critical C^ and periods for the models are given in Tabled 
In all the models we used a normalised 77 = 0.5, Cq 
and ilo for the thin and thick convection zone models are 
given in Table |21 and the parameters for the model used 
in Fig. Hare VL^ = -0.3, ri„ = 0.2 and = 0.3 Ca = 8. 

3. From the model to observations 

For converting the possible spot pattern from the model 
calculation into synthetic photometric observations, we 
first have to decide which value of the magnetic pressure 
on the stellar surface results in a spot. For doing this a 
three temperature model was chosen in which the values of 
magnetic pressure that are > 70 % of the maximum value 
are considered to form the "umbra" and the values < 70 % 
and > 30 % of the maximum form the "penumbra" . The 
values < 30 % of the maximum denote the unspotted sur- 
face. For investigating the long-term changes in the spot 
strength the maximum value of the magnetic pressure was 
taken from the whole run, not from the individual maps. 

A typical star showing flip-flops is a cool giant or a 
zero age main sequence object. For describing the realis- 
tic spot temperatures on such stars, 5000 K was chosen as 
the unspotted surface temperature and 3500 K & 4250 K 
as the umbral and penumbral temperatures, respectively. 
After the assignment of the spot temperatures to the mag- 
netic pressure maps, synthetic light-curves were calculated 
from the maps. The limb-darkening coefficient from Al- 
Naimyi H1978|l for 5000 K at the central wavelength of the 
Johnson V band was used for all the three temperatures. 
Fig. n shows the magnetic pressure map obtained from 
the dynamo calculations, the corresponding spot config- 
uration and the synthetic light-curve calculated from the 
temperature map. 

4. Results 

The results from the thick and thin convection zone 
flip-flop models, that were first discussed in Elstner & 
Korhonen H2005fl . have been converted into light-curves 
as described in the previous section. The model parame- 
ters are given in Table |3 

In Fig. 121 an example of a sequence of temperature 
maps exhibiting a migrating spot pattern and a flip-flop 
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a) Magnetic pressure b) Temperature ^ c) Photometry 




Fig. 1. From the model calculations to the photometric 
observations, a) The magnetic pressure on the surface ob- 
tained from the model calculations. The brightest areas 
correspond to the largest magnetic pressure, b) The tem- 
perature map showing the spot configuration correspond- 
ing to the magnetic pressure distribution. In the map the 
black area has a temperature of 3500 K, grey corresponds 
to the temperature of 4250 K and white is the unspot- 
ted surface with the temperature of 5000 K. c) The nor- 
malised synthetic light-curve corresponding to the tem- 
perature map. 

Table 2. The parameters for the models investigated in 
this paper. Table gives the name used for the model in the 
text, location of the inner boundary, the strength of the 
differential rotation (for solar differential rotation flo = 1), 
dynamo- number Ca (Eq.|31), the energy density in equipar- 
tition units EO for the component m = and El for m — 1, 
the period of the oscillation in diffusion times and the mi- 
gration period, also in diffusion times. 



name 


Tin 






EO 


El 


PO 


PI 


thick 


0.4 


0.12 


10 


0.06 


0.06 


0.3 


0.2 


thin 


0.7 


0.11 


21 


0.1 


0.4 


0.23 


0.46 



event is shown. The features are symmetric with respect to 
the equator. It is clearly seen that in this thick convection 
zone model the spot shift in the flip-flop is ^ 90° and not 
180°, as is more commonly seen in the observations. 

For investigating the long-term photometric behaviour 
obtained from the models, calculations were done starting 
around 80 diffusion times, running 50000 timesteps (one 
timestep is 1.9 • 10~^ diffusion times or approximately 1/5 
days for our models) and calculating a map of magnetic 
pressure at the surface every 100 steps. These maps were 
then converted into temperature maps and synthetic light- 
curves, and the light-curves were plotted against time to 
see the long-term behaviour. Figs.lSt^ & b show the calcu- 
lated light-curve behaviour for three different inclination 
angles for the thick and thin models, respectively. 

In the case of the thick model (Fig. we see that 
the minimum magnitude is very strongly modulated dur- 
ing the flip-flop cycle, whereas the maximum magnitude 
is much less affected, except when viewed from small in- 
clinations. The minimum and maximum in the maximum 
magnitude occur around the time of the minimum and 
maximum in the minimum magnitude. This behaviour is 
the same as seen by Fluri & Berdyugina (2004") in their 
first case (sign-change of the axisymmetric component). 
The inclination of the star affects mainly the amplitude of 
the variation, as it determines how close the spots are to 
the centre of the visible disk (location of the maximum ef- 



J 6 11 




Fig. 2. Temperature maps in a cylindrical projection 
showing a migrating spot pattern and a flip-flop event 
for the thick model. The colour coding is the same as in 
Fig. Each maps has 500 time steps in between them. 
Note that the small spot in the lower right corner is for 
enabling the correct plotting of the colours. 

feet on the light-curve). In the case of high-latitude spots, 
as seen in these models, this effect totally dominates over 
the other inclination effect, i.e. how much of the spots on 
the "southern" hemisphere are visible. The changes in the 
inclination affect the behaviour of the maximum magni- 
tude more strongly than the minimum magnitude. 

The behaviour seen in the thin convection zone model 
(Fig. Oh) is different from the one seen in the thick model. 
The amplitude of the variation is much smaller because 
the active longitudes are further apart than in the thick 
case. It is also worth noting, that in this case the variation 
in the maximum and minimum magnitudes is very differ- 
ent from the thick case; here the maximum of the max- 
imum magnitude occurs near the minimum of the mini- 
mum magnitude. There is a small shift towards the later 
timesteps for the maximum of the maximum magnitude in 
comparison to the minimum of the minimum magnitude. 

5. Discussion 

5.1. Dependence of the solution on the model 
parameters 

The differential rotation needed for a flip-flop solution 
with a period of about 5 years is rather small. For a solar 
sized star the latitudinal difference in (Afi) should be 
about 10% of the solar value, independent of the global 
rotation. For i?, = 2Rq the Ail can be about 40%. The 
situation seems similar in all cases, i.e. for isotropic a ten- 
sor used together with a differential rotation depending on 
the distance to the rotation axis and with an anisotropic 
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a) 



inclination =35 



nclination =60 



b) 





inclination =85 



10000 20000 30000 40000 50000 
Timestep 

inclinotion =35 





10000 20000 30000 40000 50000 
Timestep 



10000 20000 30000 40000 50000 
Timestep 




10000 20000 30000 40000 50000 
Timestep 



nclinotion =60 



10000 20000 30000 40000 50000 
Timestep 

inclinotion =85 



10000 20000 30000 40000 50000 
Timestep 



Fig. 3. The normalised magnitudes, calculated from the thick (a) and thin (b) convection zone flip-flop models, plotted 
against the timestep in the model for three different inclination angles. The series consists of synthetic light-curves 
calculated from 500 surface maps with 100 timesteps in between each map. The length of the whole series is thus 
slightly longer than 10 000 days (~28 years). The starting time for the series is around 80 convective turn-over times. 
The magnitudes have been normalised to the mean magnitude of the whole set. 
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Fig. 4. Magnetic energy in axisymmetric (m=0) and non- 
axisymmetric (m=l) components for a model with Vc = 
0.3 and isotropic a. The pole to equator difference of 
is 30% of the solar value. The mixed component flip-flop 
solution survives only 5 diffusion times. 



a tensor used with both a solar-like rotation law and axis 
distance dependent rotation law. An example of the time 
evolution of the magnetic field energy in the m=0,l com- 
ponents with isotropic a, AO 30% of the solar value, and 
a rotation law depending on the axis distance (cf. Moss 
2005) is shown in Fig. ^ As can be seen, with relatively 
strong differential rotation the non-axisymmetric compo- 
nent, m=l, is initially excited, but the mixed component 
solution survives only approximately 5 diffusion times. 

Increasing the diffusivity gives a strong flip-flop phe- 
nomenon also for higher pole to equator differences of 
Q, but with a smaller flip-flop period. The flip-flop so- 
lutions appear preferential for a positive a in the north- 
ern hemisphere in the case of radially increasing Q at 
the equator. This leads to a poleward migration of the 
spots. Observations indicate solar-like equatorward mi- 
gration pattern in solar-like stars (see e.g. Katsova et 
al. 12003(1 . but there is also a detection of poleward mi- 
gration of the spots in the RS CVn binary HR 1099 (Vogt 
et al.HnMI Strassmeier & Bartus i2000)l . 



5.2. Comparing the synthetic light-curves with the 
observations 

Many active stars show long-term light-curves where time 
periods with small and large amplitude in the photometry 
alternate, as seen in our flip-flop models. For the behaviour 
seen in the thick model (Fig.|3^) a good stellar counterpart 
is DX Leo (see e.g. Messina & Guinan l!j0()2ll . It is easier to 
find stellar counterparts for the thin case (Fig. |3ti). Some 
examples, like LQ Hya and EI Eri, can be seen for instance 
in Olah & Strassmeier (;20^. The fact that the thin case 
seems to be dominant implies that the flip-flops where the 
spots shift 180° are more common. 

As seen in the case of FK Com (Olah et al l2005|l . some 
stars can show both 90° and 180° shifts in the spots dur- 
ing a flip-flop event. For investigating what alternating 
90° and 180° flip-flops would look like in the long-term 
photometry, we combine the calculated light-curves from 
the models showing 90° (from the thick model) and 180° 
(from the thin model) flip-flops, taking alternatively one 
90° flip-flop and one 180° flip-flop. The result of combining 
the two types of flip-flops is shown in Fig. Et. The long- 
term light-curve behaviour produced by this is similar, but 
not identical, to the second case of Fluri & Berdyugina 
1)2004(1 . which shows alternatively small and large ampli- 
tude changes that are symmetric with respect to the mean 
magnitude, i.e. the maximum of the maximum magni- 
tude occurs at the time of the minimum of the minimum 
magnitude. Our combination of the models with different 
spot shifts in the flip-flop event is symmetrical only during 
the smaller amplitude phase. During the larger amplitude 
phase the minimum of the maximum magnitude occurs 
at the time of the minimum of the minimum magnitude. 
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Fig. 5. a) The combination of the synthetic light-curves from the 90° and 180° flip-flops running approximately 9000 
days, every other flip-flop is from different model. The 180° flip-flop calculations were moved slightly up for a more 
continuous plot, b) Light-curve of a Gem (from Henry et al. I.1995') shows indications of similar behaviour as seen in 
the panel a. Note the longer time span of the model (a) in comparison to the observations (b) . 



For example, the light-curve of a Gem (see Fig.|3)) shows 
indications of this kind of behaviour. 

It is often difficult to see in the real observations the 
activity pattern caused by the flip-flop behaviour. This is 
partly due to the years long time series of observations 
needed to see the pattern and partly due to the solar-like 
cyclic changes in the over-all activity level of many active 
stars. Quite drastic changes in the brightness of some stars 
can be seen on top of the possible flip-flop signature (see 
e.g. HK Lac and IL Hya in Olah & Strassmeier I2002|l and 
these large changes are likely to mask the patterns caused 
by the flip-flops. 

6. Summary 

We have synthesised photometry from Dynamo calcula- 
tions exhibiting flip-flop behaviour. This was done for in- 
vestigating the long-term changes in the photometric be- 
haviour seen over several flip-flop cycles. On the whole, 
the activity patterns discussed in this paper imply flip- 
flop phenomenon and stars showing these patterns should 
be further investigated for checking if they really show flip- 
flops. A statistically significant sample of stars is needed 
for deeper understanding of this phenomenon. Further, 
more effort should be put to measuring the Ail, merid- 
ional flow, and latitude migration of the spots on different 
types of stars. All these parameters have important impli- 
cations for the dynamo calculations. 
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